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Abstract 

We have observed /3 + -delayed a and pa emission from the proton-rich nucleus 21 Mg produced at the ISOLDE 
facility at CERN. The assignments were cross-checked with a time distribution analysis. This is the third 
identified case of (3pa emission. We discuss the systematic of beta-delayed particle emission decays, show 
that our observed decays fit naturally into the existing pattern, and argue that the patterns are to a large 
extent caused by odd-even effects. 
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1. Introduction 

Beta-delayed particle emission is an important 
decay mode for exotic nuclei and allows many as¬ 
pects of nuclear structure to be probed, see the two 
recent reviews mm for a comprehensive overview. 
We report here the first observation of (3a emission 
as well as the rare /3pa emission from the nucleus 
21 Mg. Based on these observations we have iden¬ 
tified systematic patterns in the occurence of beta- 
delayed particle decays in proton-rich nuclei. We 
shall present and discuss these as well. 

A detailed description of beta-delayed particle 
emission must include consideration of local nuclear 
structure effects, but its occurence is in general 
dominated by the available energy, i.e. the differ¬ 
ence between the Q^-value and the particle separa¬ 
tion energy. As is well known, for an isobaric chain 
with mass number A the Qp values will increase and 
the proton and neutron separation energies decrease 
as one moves from the beta stability line towards 
the driplines (modulated for even A by the pair¬ 
ing term). The a particle separation energy tends 
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for light nuclei to be minimal for N = Z nuclei, 
but the minimum moves towards more proton rich 
nuclei and reaches the proton dripline at A ~ 50. 
This causes a clear pattern for beta-delayed multi¬ 
particle emission, with /32p and /33p taking place 
close to the proton dripline, /32n, /33n etc starting 
from about halfway to the neutron dripline, while 
/32a is seen from A = 8, 9,12 nuclei close to stabil¬ 
ity. (To the extent that these decays are sequen¬ 
tial one can of course regard them as /3a decays 
to the unstable A = 5, 8 nuclei.) Similar patterns 
appear in beta-delayed single-particle emission al¬ 
though exceptions occur for the very light nuclei 
such as the large P n values for N = 10 nuclei and 
the (3a emission from neutron-rich N-isotopes. 

We focus first on the multi-particle (3pa decay 
and return in the discussion to the general patterns 
of beta-delayed particle emission. 

2. Experimental results 

The 21 Mg activity was produced at the ISOLDE 
facility at CERN by a 1.4 GeV proton beam im¬ 
pinging upon a SiC target. The produced atoms 
were extracted, laser ionized, accelerated to 60 keV, 
led through a mass separator into the experimen- 
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Figure 1: (Color online) A E-E plot from the Gas-Si tele¬ 
scope with the gas channel on the vertical axis and the 
deposited energy in the silicon detector on the horizontal 
axis. The scaled stopping powers for a’s and protons are 
shown on top of the data in solid and dashed red, respec¬ 
tively. The graphical cut used for the a-particles is shown 
with the dashed black closed line and the cut for the protons 
is shown by the solid black closed line. The events with 3.18 
MeV in the silicon detector and high energy deposition in 
the gas are due to a contamination of 148 Gd. 

tal set-up, and implanted in the window of a gas-Si 
telescope opposed by a Si(DSSSD)-Si telescope. A 
full account of the experimental procedure is given 
in 0- The collected source also contained a sub¬ 
stantial amount of 21 Na. 

The data from the Gas-Si charged particle tele¬ 
scope are presented as a A E-E spectrum in Fig. 
[l] Rescaled stopping powers [1] for a particles and 
protons (evaluated for silicon, but representing the 
total energy loss in the collection foil, the gas detec¬ 
tor and the Si dead layer) are drawn in the figure 
and match the data well, indicating the presence 
of /3a decays on top of the previously established 
[SI [S] /3p. The /3-particle component in the lower 
left corner of Fig. [l] overlaps with protons below 
1150 keV and a-particles below 700 keV making 
particle identification difficult at low energy. The 
a-particles are stopped in the DSSSD and cannot 
be separated there from the more intense proton 
branches. 

The a-particle spectrum extracted from the gas- 
telescope by applying the gate drawn as a dashed 
black closed line in Fig. [T| is shown in Fig. [2] Apart 
from a remaining background component at low en¬ 
ergy five a branches can be identified in the spec¬ 
trum. The a 2 , 03 , 04 and as lines naturally fit 
into the 21 Mg scheme put forward in [3j as /3- 
delayed single a branches (see [5j for the full decay 


Figure 2: (Color online) Singles a spectrum extracted from 
the Gas-Si telescope. The inset in the top right corner shows 
the DSSSD proton spectrum in coincidence with the low en¬ 
ergy ai line. 

scheme). The line, with measured laboratory 
energy 714(12) keV, does not fit with a transition 
between known levels in 21 Na and 17 F. However, it 
does agree with a known a-particle transition from 
20 Ne to 16 0 observed in the decay of 20 Na [7] with 
a laboratory energy of 714(4) keV. 

A conclusive particle identification for ai was not 
possible from the A E-E plot, but strong support 
for the above assignment comes from the observa¬ 
tion of a coincident line in the DSSSD detector, 
assigned to be the preceding proton. This proton 
branch P 3 (the numbering is chosen to be consistent 
with the full data set discussed in m is displayed 
as the inset in Fig. [2j From the measured energy 
we deduce -E C m(P 3 ) = 919(18) keV which leads to 
the interpretation of ai and P 3 as being due to /3pa 
decay of 21 Mg through the 5/2+ isobaric analogue 
state (IAS) at 8.975 MeV in 21 Na via proton emis¬ 
sion to the 5.621 MeV 3“ level in 20 Ne and finally 
a emission to the ground state of 16 O. The total 
branching ratio of this decay branch is found to be 
1.6(3) • 10 -4 . This proton branch from the IAS has 
not been observed earlier and a-emission from ex¬ 
cited states of 21 Na have only been reported in one 
earlier experiment [3]. 

2.1. Time distribution analysis 

As mentioned above our data are contaminated 
by 21 Na, other small contaminants could in prin¬ 
ciple also be present. The observed /3a and /3pa 
branches are quite weak, so a cross-check of the as¬ 
signment is valuable. This is done by considering 
the time distribution of the events. 
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Several factors influence the time distribution of 
the recorded 21 Mg events, see |3] for an exhaustive 
discussion. We shall use as reference the experimen¬ 
tal time distribution recorded for events within the 
proton gate in Figure[l]and with energy above 1150 
keV. The energy gate ensures that the reference dis¬ 
tribution only contains protons from the decay of 
21 Mg. The halflives of 21 Mg and 21 Na, 122(2) ms 
and 22.49(4) s [9], differ greatly as do the corre¬ 
sponding time distributions. Other contaminants 
are also expected to differ from 21 Mg. 

Some of the (3a branches have quite low statis¬ 
tics and we therefore compare their time distribu¬ 
tion directly to the reference distribution. This can 
be done very efficiently with the empirical distribu¬ 
tion function (EDF) statistics pTO] that give pow¬ 
erful goodness-of-fit tests. The basic principle is to 
compare the shape of the data sample to the refer¬ 
ence shape by measuring the distance between the 
two cumulated distributions. For experimental and 
reference distributions with values EDFi and Fi in 
bin i, the three most frequently used EDF statistics 
are m Kolmogorov-Smirnov 

D = \fN max |EDFi — Fi\, 

i 

Cramer-Von Mises 

W * 2 = Nj2(ZDF t - F x ) 2 Pl 

i 

and Anderson-Darling 

3 2 _ at (EDFi - Fi) 2 pi 

A N ^ fi(l-fi) ’ 

where N is the total number of counts and pi is 
the probability to be in bin i in the reference dis¬ 
tribution. The second column of Table [I] gives the 
95% confidence levels for the three EDF statistics 
obtained through Monte Carlo simulation, values 
below these levels indicate the time distribution for 
the different lines are consistent with the one of 
21 Mg. More details on the confidence levels are 
given in [3j. 

The EDF test results in Table [T] show that all 
lines, except for on, agree with the reference distri¬ 
bution. The agreement is particularly good for the 
strongest line, a 3 . As mentioned above there is a 
contamination of /3-particles in a\ that come from 
both 21 Na and 21 Mg. We would therefore expect 
the time distribution for a\ to be mainly that of 
21 Mg with a small component of 21 Na. The EDF 
tests are sufficiently sensitive to see the effect of 


Table 1: EDF goodness-of-fit tests of the time distribution 
of the individual observed lines. The first column denotes 
the test, the second column gives the 95% confidence level 
(obtained through Monte Carlo simulations) for having the 
21 Mg time distribution. 
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1.08 
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0.33 
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2.49 

1.46 

3.61 

0.64 
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the small 21 Na contribution. We expect, and do 

observe, that the upper part of the a\ distribution 
has smaller contamination level. The fact that the 
coincident p 3 distribution is consistent with being 
from 21 Mg implies that we can safely assign the 
/3pa transition, as well as all /3a transitions, to the 
decay of 21 Mg. 


3. Discussion 

3.1. Other /3 pa cases 

The /3p a decay mode is very rare as described 

in [TJ with only two previously established cases: 

9 C and 17 * * * * * * * Ne. For two further candidates, 13 0 and 

23 Si, the decay mode has not been seen so far. Most 

searches have concentrated on seeing particle emis¬ 

sion from the IAS in the beta-daughter due to the 

large beta-strength to this state. 

The case of 9 C is special in that all states popu¬ 

lated in the beta-daughter 9 B break-up into two 

a-particles and a proton, see [12] and references 

therein. This could be presented as a 100% branch¬ 
ing ratio for /3pa or /3ap decay to 4 He, but the 
decays of the A = 9 nuclei are special in several as¬ 
pects HI El and are not typical for this decay mode. 

Although /3pa has not been observed so far for 
13 O it must occur since /3-decays to the IAS in 13 N 
have been observed m and close to half of the 
IAS decays are known from reaction experiments 
m to go via proton-emission to a unbound states 
in 12 C or a-emission to proton unbound states in 
9 B. Actually, the final state in both cases will be a 
proton and three a-particles which makes the decay 
more challenging to observe. The total branching 
ratio for the decay mode can be estimated to be 
0.9(3) • KT 4 . 

For 17 Ne both decay orderings, /3pa and /3ap, 
have been observed m with a total branching ratio 
for the decay mode of 1.6(4) • 10 -4 . All observed 
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Figure 3: (Color online) The /3+-decaying isotopes of the 
elements from N to Si. Dark squares indicate stable iso¬ 
topes. The experimentally observed beta-delayed particle 
decay modes are indicated, the /3p a decay modes for 21 Mg 
(seen for the first time here) and 13 O (see section 3.1) are 
marked in blue. 


decays proceed through the IAS in 17 F and go to 
the final nucleus 12 C. 

Adding now our observation of 21 Mg(/?pa) 16 0 it 
is striking that all cases go through an a-conjugate 
nucleus, namely 8 Be, 12 C, 16 0 and 20 Ne respec¬ 
tively. Before drawing any firm structure conclu¬ 
sions we shall consider the broader systematics of 
beta-delayed particle emission in Z > N nuclei. 

3.2. Systematics of beta-delayed decays 

Similar patterns also appear in other beta- 
delayed particle decays (see mm for more data and 
for references to the original work). One closely re¬ 
lated example is /3a decays that occur for all bound 
A = 4n, T z = — 1 nuclei up to A = 40: 8 B, 12 N, 
20 Na, 24 Al, 28 P, 32 C1, 36 K and 40 Sc. The /3p de¬ 
cays are well established I] to occur strongly in 
A = An + 1, T z = —3/2 nuclei. The /32p decays of 
22 Al and 26 P and the /33p decay of 31 Ar [1B| also all 
end up in an a-conjugate nucleus. The decays ob¬ 
served for the elements N to Si are shown in figure 
[3] Note that the /3a and /3p modes are not marked 
explicitly when the /3pa or /32p modes also occur. 

In the following we shall argue that the observed 
patterns are likely (except for the very lightest nu¬ 
clei) to be related to odd-even effects rather than 
a-cluster structure. We start by considering the 



Figure 4: (Color online) Qec values for nuclei up to Ge. 
Filled symbols indicate nuclei with even Z, open symbols 
odd Z. Squares indicate nuclei with even N, circles odd N. 
The isospin values given are for the even-even nuclei. The 
green and red symbols correspond to T = 1 ( T z = 1) and 
T = 4. The dashed (dotted) line is a liquid drop estimate 
°f Qec for nuclei with T z = —1/2 (—3/2). The full lines 
indicate quartets of Al-Si nuclei. 


systematics of Qsc-values for nuclei with Z > N 
as illustrated in figure |4j Even though many effects 
contribute to the masses in this region, a liquid drop 
estimate reproduces the trend of Qec for the odd- 
A nuclei with T z = —1/2 (dashed line) where only 
the Coulomb term enters, as well as for T z = —3/2 
(dotted line) where the asymmetry term also con¬ 
tributes. Note that Qec in the latter case varies 
little for A between 25 and 50. 

The experimental data show that the Q ^c-values 
are roughly the same for each “quartet” of four nu¬ 
clei that, as illustrated in the left panel of figure [5j 
have proton and neutron numbers (Z,N), (Z,N+1), 
(Z-1,N), (Z-1,N+1) where both Z and N are even. 
This is pronounced for quartets where the even-even 
nucleus has T = 1, and holds to a lesser degree also 
for T = 2 for mass numbers up to 40. The reason 
for this is that the two odd-A nuclei are at the same 
distance from the beta-stability line and therefore 
have about the same Q-value, as also shown by the 
liquid drop estimate. Without a pairing term in the 
liquid drop formula the Q-value for the even-even 
nucleus would be larger and the odd-odd smaller, 
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Figure 5: Left part: The quartet of nuclei with similar Q- 
value. Right part: the favoured beta-delayed decay modes. 
See the text for details. 


but the odd-even effects counteracts this and as can 
be seen from figure [4] the magnitudes are even re¬ 
versed for most nuclei. For the quartet with T = 1 
(and T z = —1) the odd-odd nucleus has N = Z 
and is therefore extra bound, this happens to re¬ 
sult in Q-values that are almost the same for all 
four nuclei. The quartets are indicated in figure [3] 
by thicker lines. 

The observed decay patterns now follow from the 
energetics and are illustrated in the right panel 
of figure [5] The /3a decays should occur in odd- 
odd nuclei, since they have slightly higher Q -values 
and the daughter alpha particle separation energies 
tend to be smallest here. The /3p decays need low 
proton separation energies in the daughter nucleus 
and therefore are more prominent for even Z, start¬ 
ing (as one goes from stability towards the proton 
dripline) in an even-odd nucleus. The /3pa decay 
should be favoured in even-odd nuclei, and /32p and 
/33p decays should occur in odd-odd and even-odd 
nuclei, respectively, by extending these arguments. 

Experimentally, the /3a, /3p and /3pa decays ap¬ 
pear first in the quartets where the even-even nu¬ 
cleus has T = 2 and the odd-odd nucleus T = 1, 
but /3p occurs also in 59 Zn, 65 Ge and heavier nu¬ 
clei. The beta-delayed multi-proton decays appear 
in more exotic nuclei, but it is noteworthy that /3a 
in these nuclei only has been observed in the odd- 
odd 22 Al. 

Similar patterns can be expected for /3 _ -delayed 
particle decays, although the grouping of Q-values 
is less pronounced here. The /3a and /3na decay 
modes will in general occur further away from the 
beta-stability line. 


Our study of the decay of 21 Mg has given the first 
evidence for the occurence of the /3a and /3pa de¬ 
cay modes in this nucleus. The assignment of these 
decay modes to 21 Mg has been verified through sta¬ 
tistical tests of the time distribution of the events. 
The occurence of these decay modes in the /3 + decay 
of 21 Mg fits naturally into the systematics of previ¬ 
ously observed /? + -delayed decays. We presented a 
brief overview of this systematics and argued that 
it can be explained by the variation in decay energy 
due to odd-even effects and that there is no need 
to invoke specific structure effects such as alpha¬ 
clustering in spite of a-conjugate nuclei occuring 
often as final state nuclei. 

This interpretation can be tested when new in¬ 
stances of these exotic decays are discovered. The 
/3pa decay mode may not occur in heavier XI = 
—3/2 nuclei than 21 Mg (the Q-value becomes more 
than 10 MeV in 61 Ge, but the Coulomb barrier for 
a-particle emission is substantial then), but may be 
found also in the T z = —5/2 nuclei 23 Si, 27 S, 31 Ar 
etc. If found in 20 Mg it may help to quantify the 
15 0(a,7) 19 Ne reaction rate [17]■ A general overview 
of which energetically allowed decays have not yet 
been observed was given in |T[. 
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